Spin polarized transport through a single-molecule magnet: 
current-induced magnetic switching 
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Magnetic switching of a single-molecule magnet (SMM) due to spin-polarized current is investi- 
gated theoretically. The charge transfer between the electrodes takes place via the lowest unoccupied 
molecular orbital (LUMO) of the SMM. Generally, the double occupancy of the LUMO level, and 
hence a finite on-site Coulomb repulsion, is allowed. Owing to the exchange interaction between 
electrons in the LUMO level and the SMM's spin, the latter can be reversed. The perturbation 
approach (Fermi golden rule) is applied to calculate current-voltage characteristics. The influence 
of Coulomb interactions on the switching process is also analyzed. 

PACS numbers: 72.25.-b, 75.60.Jk, 75.50.Xx 



I. INTRODUCTION 

Single-molecule magnets (SMMs) are molecules with 
a relatively large net spin moment (corresponding to 
the spin number S) and a significant uniaxial magnetic 
anisotropy.^ As a result, behavior of SMMs resembles 
much that of superparamagnets, and at low tempera- 
tures the molecules become trapped in one of the two 
metastable spin states j ± S) Owing to this bistabil- 

ity, SMMs seem to be a suitable base for memory cells of 
future information storage and processing technology. 7,8 
Apart from this, SMMs can possibly become basic com- 
ponents of the molecular-based spintronic devices* 9 - 

Electronic transport through an individual SMM 
has been demonstrated experimentally only very re- 
cently ) 11 : 12 : 13 attracting also a significant theoretical at- 
tentioni 14 i 15 i 16 i 17 i 18 An important issue in this context is 
the question of how the molecule's spin can be switched 
between the two stable states by means of spin-polarized 
current. This question is important not only from the 
purely fundamental reasons, but also from the point of 
view of possible applications of SMMs in various magne- 
toelectronic devices, and particularly as memory cells. 

It is already well known that when a spin-polarized 
current flows through a magnet, some amount of spin 
momentum can be transferred from the current to mag- 
netic body, 19 leading effectively to a spin-transfer torque. 
This additional torque may lead to magnetic switch- 
ing or current-induced precessional states. It has been 
shown recently, that exchange interaction between spin- 
polarized current and a SMM embedded in the bar- 
rier of a magnetic tunnel junction can lead to rever- 
sal of the molecule's spin^ The model considered there 
was simplified as the current was not flowing through 
the molecule, but rather directly between magnetic elec- 
trodes. However, the tunneling electrons could inter- 
act with the SMM via the exchange coupling, leading 
to switching of the SMM. The main objective of the 
present paper is to investigate theoretically a more re- 
alistic mechanism of SMM's switching, when the spin- 
polarized current flows directly through the molecule 
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FIG. 1: (Color online) Schematic representation of the system 
under consideration in the nonequlibrium state, i.e. when a 
finite bias voltage V is applied, eV = fj,L — fJ,R, where Hl{^r) 
denotes the electrochemical potential of the left (right) elec- 
trode. Two collinear magnetic configurations of the leads' 
magnetic moments: parallel (black solid arrows) and antipar- 
allel (grey solid arrows) are also indicated. 



(molecular single-electron transistor geometry). In the 
model assumed, the current flows via the lowest unoccu- 
pied molecular orbital (LUMO) of the SMM. We restrict, 
however, our consideration to the case of the sequential 
transport regime. The results clearly show that trans- 
port of electrons through the LUMO level can lead to 
magnetic switching of the molecule, when the electrons 
in the LUMO level interact via exchange coupling with 
the spin moment of the inner core of the SMM. 

In Sec. II we present the model and basis of the theoret- 
ical analysis of transport characteristics. Numerical re- 
sults are presented in Sec. III. These results clearly show 
the possibility of magnetic switching induced by current 
pulse. Summary and final conclusions are in Sec. IV. 



II. MODEL AND THEORETICAL 
DESCRIPTION 

The model to be considered, see Fig. [TJ consists of 
a SMM weakly coupled to two ferromagnetic electrodes 
(also called here leads). We assume that the electronic 
transport through the molecule takes place only via the 
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LUMO level of energy s (measured from the Fermi level 
of the leads at equilibrium). This level is assumed to 
be exchange-coupled to the SMM's spin. 8 Only collinear, 
i.e. parallel and antiparallel configurations of the leads' 
magnetic moments are considered, and these magnetic 
moments are assumed additionally to be parallel to the 
magnetic easy-axis of the molecule (along the axis z in 
Fig. p. 

The complete Hamiltonian of the system may be writ- 
ten as H = Hsmm + H e i + Ht- The first term describes 
the SMM and is assumed in the form, 

Hsmm = -(d+ AD 1 cic a + AD 2 c\c 1 c\c^S 2 
<7={T,1} 

+ e4 c <7 + Uc \ c \ c \ c l 

~2 E °W sr 4 £ V> (!) 

where S is the SMM's spin operator, er is the Pauli 
spin operator for electrons in the LUMO level, and 
cJ.(c CT ) is the relevant electron creation (anihilation) op- 
erator. Apart from this, U denotes the Coulomb energy 
of two electrons of opposite spins in the LUMO level, 
whereas J is the exchange coupling parameter between 
the SMM's spin and electrons in the LUMO level. The 
parameter D is the uniaxial anisotropy constant of a 
free molecule, while AD\ and AD 2 describe corrections 
to the anisotropy when the LUMO level is occupied by 
one and two electrons, respectively^ The perpendicular 
anisotropy terms have been omitted as irrelevant for the 
effects described here. Apart from this, we neglect intrin- 
sic spin relaxation, eg. that due to spin-phonon coupling. 
The only spin relaxation taken into account is that due 
to coupling of the molecule to the electrodes. This is 
justified as spin relaxation due to electronic processes as- 
sociated with coupling of the dot to external leads is the 
dominant one. 

The next term of H describes the ferromagnetic elec- 
trodes, 
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F q n qt n q 
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where L(R) stands for the left (right) lead. The leads 
are characterized by non-interacting electrons with the 
dispersion relation g£ , where k denotes a wave vector. 
9 



In the equation above, aj. CT and a 



are the relevant an- 
nihilation and creation operators, respectively. 

The final term of the Hamiltonian H represents tun- 
neling processes between the leads and the molecule, 



"* = E E [ 

q={L,R}k,<T={U} 
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Tidal, 
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where Tl and Tr are the tunneling matrix elements be- 
tween the SMM and the left and right electrodes, respec- 
tively. These parameters are assumed to be independent 



of the wave vector and spin orientation. We point, that 
direct electron tunneling between the leads is excluded. 

It has been shown 8 that the Hamiltonian Hsmm, 
Eq. ([T]), commutes with the z component S* of the to- 
tal spin operator St = S + J2aa r <T <^^' c t c a' /2j where the 
second term of St represents the spin of an electron in 
the LUMO. As a consequence, if one treats Hsmm as the 
unperturbed part of the total Hamiltonian H, it is con- 
venient to numerate the eigenstates of Hsmm with the 
eigenvalues m of Sf and with the number of electrons in 
the LUMO level. Thus, the eigenstates of the SMM are 
given by: |0,m) = |0) o <8> \m) mo i, \2,m) = \ T|) ® \m) m oh 
|l,m) ± = A±||) ®|m + l/2) mo/ + B±|T) ®|m-l/2) mo/ 
for the intermediate states, and |1,±5 ± 1/2) = | T (I 
))o <g> I ± S) mo i for the fully polarized states. Accord- 
ing to our notation, |») (moi) denotes the spin state of 
the orbital (SMM). The coefficients A+ and B+ act here 
as effective Clebsch-Gordan coefficients which depend on 
the system's parameters, and have the form 



± _ ^ ^2A^ 
v — + 



± (2/jM - J)i 



2^A^r 



Jy/S(S+l) 



1/4 



(4) 
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2 x /Ae(m) v /2Ae(m) ± (2LK 1 ) - J)m 

where Ae(m) = y/DW(DW - J)m 2 + (J/4) 2 (25 + l) 2 
and L>W = D + ADl Additionally, we assume 
J ^ 0. The corresponding eigenenergies of the Hamilto- 
nian Hsmm are: e(0, m) = —Dm 2 , e(2, m) = 2e + U — 
(D + 2AD 1 + AD 2 )m 2 and e(l, m) ± = e + J/4 - (D + 
ADi)(m 2 + 1/4) ± Ae(m). The energy of the fully po- 
larized states |1,±5± 1/2) is e(l, ±S ± 1/2)+. 

To investigate the current-induced magnetic switching 
of the SMM, we analyze the relevant I-V characteristics. 
The total current flowing through the molecule can be 
written as I = (II — Ir)/2, where I a (a — L, R) denotes 
the current flowing from the lead a to the molecule, 



E E(* 



\ \n q ,m q )\n r ,m r ) p (a\ 
V la r \n q ,m q )- V"^ 



m r ,m n r ,n a 



Here, ^t' m ii\ nr ' mr i represents the rate of transitions be- 
tween the states \n q ,m q ) and \n r ,m r ), whereas P\ nq ,m q ) 
is the probability of finding the SMM in the state 
\n q ,m q ). We assume that current is positive when elec- 
trons flow from the left to right. For notational clarity, 
from now on we assume \n q , m q ) = \q), which also means 
that E 9 = En, E„ v 

To find current we need to determine first both the 

\a)\r) 

transition rates 7« and the probabilities P\ q \. Let us 
start with the transition rates. In the second order (Fermi 

golden rule) one finds ja^^ in the form, 



l)\ r ) 



<) l9> /(^ 



where the first term corresponds to electron transitions 
from the a-th lead to the molecule, while the second term 
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describes the charge transfer back to the lead a. Further- 
more, f(s) is the Fermi-Dirac distribution function, and 
W) = (2Tr/h)\(f\H T \i}\ 2 5(E f - Ei) is the rate of transi- 
tions from an initial state (i) to a final state (/). 



The final expression for the transition rates yd 
the form 



l?>|r> 



takes 



\Q)\r) _ 



:{\C° r \ 2 f(e(r)-e(q) 
+ |C*| 2 [l-/(efa)-e(r) -/*„)]} 



E r ° 

={!,?} 



M«J 
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where T% = 2n\T a \ 2 D^ is the LUMO level width ac- 
quired due to coupling of the level to the lead a, and 
D" denotes the spin-dependent density of states (DOS) 
at the Fermi level in the a-th electrode. These param- 
eters will be used in the following to describe strength 
of the coupling between the SMM and leads. It is con- 
venient to write as r" = T a (l ± P a ), where r Q = 
(r? + )/2, and P a is the spin polarization of the lead 
a, P a = {D% - D1)/(D% + D°L). Here a = +(-) corre- 
sponds to spin-majority (spin-minority) electrons. In the 
following we assume that the couplings are symmetric, 
r L = Ir = T/2. Finally, in Eq. © \C a qr \ 2 = \{r\cl\q)\ 2 
together with \C^ q \ 2 = \{r\c a \q)\ 2 constitute basic selec- 
tion rules for transitions between neighboring molecular 
states. The transition is allowed only when the charge 
state of the SMM is changed by one and the change in 
the total spin satisfies AS* = ±1/2. Assuming that the 
SMM is initially saturated in the state |0, — S), see Fig.[2l 
one may expect that at a sufficiently large voltage the 
molecule can be switched to the final state |0, S). The 
switching process corresponds then to the reversal of the 
SMM's spin via all the intermediate states. 

The probabilities P\ q ^ (see Eq. (JSJ)), are obtained from 
the master equations 



dP\ 
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forn 9 = 0,lW, 2 and m q E {-S- 1/2, 5+1/2) (we recall 
the definition, \n q ,m q ) = \q)). In the present paper, we 
assume that the voltage is augmented linearly in time, 
V — ct, with c denoting the speed at which the voltage 
is increased. The corresponding time scale, however, is 
much slower than that set by electronic transitions. The 
relevant boundary conditions for the probabilities P\ q ) 
are: P\ Q _ S) (V 
\1,-S). 



0) = 1 and P ]q) (V = 0) = for \q) ^ 



III. NUMERICAL RESULTS AND DISCUSSION 

The results have been computed for an octanu- 
clear iron(III) oxo-hydroxo cluster of the formula 
[Fe 8 2 (OH)i 2 (tacn) 6 ] 8+ (shortly Fe 8 ). Such a free- 
standing molecule has spin corresponding to S — 10. The 
following values of the molecule parameters are taken: 
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FIG. 2: (Color online) The energy levels of a Fes molecular 
magnet for two values of the Coulomb parameter U. The 
black bold circle indicates the initial state |0, — S), and Ae = 
0.36 meV is the energy gap corresponding to the activation 
energy for the magnetic switching. 



D = 0.292 K (D w 0.025 meV), 22 J = 0.025 meV and 
e = 0.25 meV. Since there is no clear and reliable ex- 
perimental evidence of the change in anisotropy constant 
of the Fes molecule due to extra electrons in the LUMO 
level, we assume ADi = AD2 = 0. The couplings of 
the molecule to the left and right lead are assumed to 
be the same, Tl = = 0.0015 meV. We also assume 
that both the electrodes are made of the metallic mate- 
rial characterized by the same polarization parameter P, 
P = Pl = Pr. The calculations have been performed for 
the temperature T — 0.01 K, which is below the block- 
ing temperature Tr — 0.36 K. The corresponding energy 
levels of the molecule are shown in Fig. 2. It is worth 
noting that for the parameters assumed, the ground spin 
state of the molecule attached to the leads (<Sf = ±10) is 
the same as that of a free-standing molecule (S z = ±10). 
Nevertheless, for a sufficiently low energy of the LUMO 
level, which can be controlled for instance with a gate 
voltage, the ground state of the molecule attached to the 
leads can correspond to Sf — ±21/2 (the molecule with 
one extra electron on the LUMO level). 

Fig. [3] presents the average (S^) and current / flowing 
through the system for different values of the Coulomb 
parameter U in both parallel and antiparallel magnetic 
configurations of the leads. It can be noted that the 
reversal of the SMM's spin occurs only in the antiparallel 
configuration, whereas in the parallel configuration all 
molecular spin states gradually become equally probable. 
As a consequence, {Sf) — > and the magnetic switching 
is not observed. This is a consequence of the left/right 
symmetry of the molecule's coupling to external leads - 
similarly to the absence of spin accumulation in tunneling 
through a metallic nanoparticle in the parallel magnetic 
configuration. However, such a symmetry is absent in 
the antiparallel configuration, and accordingly the spin 
states of the molecule become unequally occupied, which 
in turn results in spin reversal. 

The corresponding current-voltage characteristics, 
shown in Fig.[3]by the dashed lines, reveal features (steps 
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FIG. 3: (Color online) The average value of the total spin 
{St ) (solid line) and the current I flowing through the system 
(dashed line) in the case of parallel (P) and antiparallel (AP) 
magnetic configurations for different values of the Coulomb 
parameter U. The other parameters are: Pl = Pr = 0.5 and 
c = 1 V/s. 
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FIG. 4: (Color online) Differential conductance dl/dV for 
two collinear, i.e. parallel (P) and antiparallel (AP) magnetic 
configurations. The parameters as in Fig. [3] 



in the current) which are directly related to the reversal 
process. In fact, each step corresponds to a certain type 
of transitions between neighboring molecular spin levels. 
Consider for instance the main features of the I-V plots 
shown in Fig. [3| The plots are almost the same as long as 

V < 0.72 mV, i.e. for eV below the energy activating the 
reversal process. The first step corresponds to the transi- 
tion between the states |0, —10} and |1, —21/2) (see also 
Fig. [2]). As voltage increases further, the next step ap- 
pears due to transitions between the states |0, —10) and 
|1,— 19/2} - , and the magnetic switching begins. Aug- 
menting voltage further, one finds certain regions of bias 
voltage, where current drops with increasing bias. This 
behavior is particularly visible in the parallel configu- 
ration, see Figs. [3] (a,c). The drop in current with in- 
creasing bias is equivalent to negative differential conduc- 
tance (NDC), see also Fig. [4] The negative differential 
conductance is a consequence of the spin blockade phe- 
nomenon^, which in turn follows from the inequality of 
the transition rates to the two levels. The suppression of 
the current starts when the system becomes energetically 
allowed to transfer from the state |1, —21/2} to the state 
|2,-10). 

The differential conductance corresponding to the I- 

V curves shown in Fig. [3] is presented in Fig. [4] The 
peaks correspond to the steps from Fig. [3J The negative 
differential conductance due to spin blockade is clearly 
seen, particularly in the parallel configuration for U > 0, 
although some weak negative peaks also occur in the an- 
tiparallel case. Evolution of the differential conductance 
with the Coulomb energy U and the bias voltage V is 
shown explicitly in Figs.[H(g)-(h) for both magnetic con- 
figurations of the leads. 



IV. SUMMARY 



In this paper we have considered electronic transport 
through a single molecular magnet attached to ferromag- 
netic leads. The molecule is characterized by a spin num- 
ber S and an additional unoccupied orbital, which be- 
comes active in transport through the molecule. 

We have shown that spin polarized electrons tunnel- 
ing through the LUMO level of a SMM can revers the 
SMM's spin when the electrons in the LUMO level in- 
teract via exchange coupling with the SMM's spin. The 
reversal starts at a certain threshold voltage correspond- 
ing to the distance between the two lowest energy levels. 
It is interesting to note, that for symmetrical systems, 
the spin reversal takes place only in the antiparallel con- 
figuration. The conductance spectra also show regions 
of negative differential conductance due to spin blockade 
effect. 
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